Universal fluctuations in tropospheric radar measurements 
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Radar data collected at an experimental facility arranged on purpose suggest that the footprint 
of atmospheric turbulence might be encoded in the radar signal statistics. Radar data probability 
distributions are calculated and nicely fitted by a one parameter family of generalized Gumbel 
(GG) distributions, G a - A relation between the wind strength and the measured shape parameter 
a is obtained. Strong wind fluctuations return pronounced asymmetric leptokurtic profiles, while 
Gaussian profiles are eventually recovered as the wind fluctuations decrease. Besides stressing the 
crucial impact of air turbulence for radar applications, we also confirm the adequacy of G a statistics 
for highly correlated complex systems. 

PACS numbers: 05.20.-y - Classical statistical mechanics; 84.40.Xb - Telemetry: remote control, remote 
sensing, radar; 05.45.Tp - Time series analysis 



Radar sensors are currently employed in a large variety 
of applications, e.g. surveillance and target identification 
[H, weather observations and forecasting 0], and geo- 
physical investigations It is therefore of paramount 
importance to have a clear control on all possible sources 
of external disturbances, which could significantly alter 
the measurements response. 

In 2006, within the activities of the GALAHAD project 
[H, a large collection of ground based radar data was 
recorded, and soon realized to constitute a unique op- 
portunity for investigating the impact of external per- 
turbations on the sought response. It was in particu- 
lar observed that occasionally, but in coincidence with 
hard windy conditions, the radar interferograms were 
affected by an overall decorrelation. Motivated by this 
unexpected finding, it was hypothesized that the atmo- 
spheric turbulence could drastically influence the radar 
propagation even over short ranges, an observation which 
however rested on purely qualitative ground. No further 
experimental realizations were in fact subsequently devel- 
oped to scrutinize the process and reach an unambiguous 
proof of concept. 

At variance, on the theoretical side, pioneering stud- 
ies on electromagnetic signal propagation through a tur- 
bulent atmosphere [f| @ date back to the 70s. Punc- 
tual fluctuations of air refractive index materialize in a 
recorded imprint, which could therefore encode an in- 
direct signature of turbulence. Indeed, the embedding 
medium can be ideally segmented in neighboring volume 
cells, each contributing to the signal according to an over- 
all trend and its own statistics. The signal impinging on 
the junction between two adjacent cells is partly reflected 
and partly transmitted and the data collected at the re- 
ceiver is sensitive to the sequence of scattering events 



occurred along the propagation path. In this sense, the 
radar echo is a global quantity indirectly representing 
the specific state of the crossed medium. A compre- 
hensive interpretative framework for such phenomena is 
however still lacking, following the inherent difficulties to 
accommodate for the intimate, highly complex nature of 
radar/air interactions. In particular, the characteristic 
of air induced fluctuations remain to be fully elucidated. 

Starting from this setting, and to gain a conclusive in- 
sight into this crucial issue, we have arranged an outdoor 
experiment, realized in a controlled environment. A k u - 
band radar was mounted on a plateau looking at an arti- 
ficial target at fixed distance. The location was chosen in 
order to reproduce as close as possible an ideal situation, 
in which the radar signal fluctuations would solely have 
depended on the propagation medium. An anemometer, 
located near the target, was used to register the wind 
velocity. 

Then, the probability distribution function of the col- 
lected signal was obtained and studied, focusing on the 
shape modification as function of the air wind condition. 
This allowed us to return a statistical characterization of 
the radar fluctuations due to the atmospheric propaga- 
tion, up to momenta of arbitrarily high order. As we shall 
demonstrate, strong wind conditions are associated to 
asymmetric leptokurtic profiles, standard Gaussian dis- 
tributions being instead recovered as the wind strength 
decreases. This observation can be cast in the form of 
a phase transition, the wind strength acting as the ex- 
ternal control parameter. To establish a clear causality 
relation between the air turbulence, here encapsulated in 
the wind velocity, and the radar signal statistics consti- 
tutes the primary goal of the Letter. Equally important, 
as an additional result of our analysis, we will confirm 
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the adequacy of the so called generalized Gumbel (GG) 
distribution as the paradigmatic PDF of global quantities 
in correlated, spatially extended system Q- 

In our experiment, a portable continuous wave stepped 
frequency fc u -band radar was installed over the excava- 
tion area of a disused quarry of Pietra Serena located 
near Firenzuola (Firenze, Italy) (22[. It operated with 
range resolution of 1 meter and unambiguous range of 
600 meter. Details on the radar system can be found in 
. The observed scene was essentially a rectangular rock 
plateau with three sides ending with a cliff [23[ . The site 
was equipped with an artificial radar target, a trihedral 
corner reflector, posed at the edge of the cliff opposite to 
the radar side (see Figure [lj . The radar-target distance 
was about 60 m. An ultrasonic anemometer, working at 
4 Hz sampling rate, installed near the radar, measured 
the wind velocity fluctuations. 

Radar data were collected at a rate of about 50 Hz, 
(sampling a 6 hours long recording window) . Mainly they 
represent the amplitude of the backscattered signal from 
the corner reflector after propagation to the corner and 
back. The data to be analyzed were arranged in a time 
series A (t) of about one million points. 

To just focus on air's impact on the propagation, all 
other possible sources of disturbances need to be iden- 
tified and subsequently excluded. In the current case, 
because of the bareness of the scene, signal fluctuations 
may depend on the oscillations of the radar tripod and/or 
the corner reflector's support, and on the intrinsic instru- 
mental noise. Mechanical oscillations can be straightfor- 
wardly characterized in the frequency domain and their 
contribution to the PDF shape was found to be negligi- 
ble. The instrumental noise is thermal noise, which is 
known to be white (SNR — 60dB). The residual signal 
as departing from the expected, a priori constant output, 
is imputed to atmospheric propagation. 

Different scale processes occur within the atmosphere, 
generally referred to as synoptic-scale processes (dozen 
of days), meso-scale processes (from day to hour) and 
micro-scale processes (less than an hour). Turbulent phe- 
nomena concern the micro scale Q , corresponding to dis- 
tances shorter than 1 km. In order to remove the compo- 
nents due to meso-scale and longer scale processes from 
the series A(t), a multiplicative detrcnding was applied. 
The "detrended" series A det (t) is determined as 



A det {t) 



in 



where Aq (i) is the "trend" obtained by a moving aver- 
age filter, which in this case is equivalent to a low-pass fil- 
ter A (t) with a threshold frequency of 5.5 TO -4 Hz. Once 
detrended, only micro-scale fluctuations are expected to 
contribute to the PDF. 

Non-stationary modulation can still affect the de- 
trended series, e.g., varying Reynolds regimes. A 60 min- 
utes long temporal window was applied sequentially, in 



order to analyze finite portions of the series, that are as 
stationary as possible, while still containing a meaningful 
number of samples for statistical analysis. Each selected 
window consists of about 10 5 points, which renders of 
some reliability the study of rare fluctuations. Different 
windows overlapped partially, for about one third of their 
length. 

Turbulent wind conditions associated to the different 
analysis windows, is quantified via the turbulent kinetic 
energy (TKE) indicator K t This latter is calculated 
from the fluctuations of the three wind velocity compo- 
nents v' x , v' y and v' z , revealed by the anemometer and 
reads: 



(2) 



where the symbol (•) stands for a time average over the 
inspected window. 

To obtain an exhaustive representation of the signal 
fluctuations, beyond the second momentum, we recon- 
structed the full PDF, within each selected time domain. 
Imagine to assume an ideal setting where interference of 
the radar signal with the wind can be neglected. Then, 
having eliminated other possible sources of systematic 
errors, one expects to recover the theoretically predicted 
PDF for a pointwise scatterer (the corner reflector) suf- 
fering from thermal noise. This is the well-known Rice 
distribution [l(| , a right skewed continuous curve, which 
approaches the Gaussian for high signal-to-noise ratio. 
Surprisingly, and at variance with what is customarily 
believed, a quite different distribution appears from our 
data analysis: the PDF of A dc t is an asymmetric, left 
skewed, distribution characterized by an exponential tail 
on the one side and a rapid falloff on the other. Such a 
significant discrepancy can be ultimately imputed to re- 
flect the non trivial interplay between propagating signal 
and turbulent dynamics of the embedding atmospheric 
medium, an effect which is not accommodated for in the 
theoretical scenario yielding to the Rice profile. 

Motivated by this working ansatz, we set down to an- 
alyze more closely the experimental PDFs. To this end, 
and aiming at a quantitative representation of the data, 
we considered a non-Gaussian, one parameter family of 
distributions, recently invoked to describe the fluctua- 
tions of global quantities in turbulent context [llj . These 



are the so-called generalized Gumbel (GG) distributions, 
orig inating from the study of extreme values statistics 
jl2| fl3j ] and more recently proposed to universally de- 
scribe the fluctuation of global quantities in correlated 
many degrees of freedom systems Q Q- I n this re- 
spect the GG constitutes the natural extension of the 
Gaussian distribution which instead apply to uncorre- 
cted variables, as a consequence of the central limit the- 
orem. In the last years, GG distributions proved to ade- 
quately adapt to distinct realms of applications, ranging 
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from e.g. the magnetization of the two dimensional XY 
model close to the Koterlitz-Thouless transition 16 1 
171 [24[ , to the statistics of the level of the Danube river 



18j . passing through the power injected in a turbulent 
flow and in electroconvetion [ijj]. Complexity and high 
correlation should be the common feature of these sys- 
tems, ending up with a GG distribution. 

The GG distributions family of curves G a is defined 
as: 



G a (x) 



t^t exp i a t 6 ( a ) ( a 



a {a)) 



a b(a)(x-s(a))~\ 



a e K+, b{a) 



d 2 
da - 



lnr(a), 



(3) 

When the real positive parameter a varies from large 
to small values, G a (x) changes its shape from that of a 
Gaussian distribution towards a more and more skewed 
PDF. This feature of the distribution family G a rendered 
it very popular in complex dynamics, when the statis- 
tics of global quantities are dealt with. The parame- 
ter a has been proposed to be inversely related to the 
correlation length of the system [l8j], and it has been 
used for measuring the number of its effective degrees of 
freedom j2(j. As a side observation, it is worth noting 
here that the GG distribution with a = ? approaches 
the well known Bramwell-Holdsworth-Pinton distribu- 
tion that was found for the fluctuations of turbulent 
power consumption 
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Back to the data collected in our experiment, properly 
rescalcd fluctuations are calculated and the normalized 
experimental histograms reconstructed. These latter are 
then interpolated via G a (x), a being adjusted as the best 
fit parameter. In Figure [2] three representative curves 
are displayed which testified on the adequacy of the pro- 
posed ansatz. The general agreement between data and 
GG distributions is remarkable. It's worth noting that 
we didn't find any correspondence between the statisti- 
cal behavior of the radar signal and TKE, approaching 
the distribution of TKE a Rayleigh distribution on every 
analyzed window. 

Even more interestingly, the values of a as obtained 
from the fitting procedure appear to be strongly cor- 
related to the quantity K t , representing the turbulent 
wind intensity within the crossed medium, for every cho- 
sen (statistically stationary) data segment. This find- 
ing, reported in the main panel of Figure [2j materializes 
in a genuine transition for a" 1 vs. K t , showing the di- 
rect footprint of turbulence in the medium on the non- 
Gaussian nature of the radar fluctuations. 

The transition between Gaussian and a non-Gaussian 
regime is rather sharp and takes place for a specific value 
of K t = K£ (~ 6 m 2 s -2 , in this case). The robustness of 



the scenario here depicted is also tested versus modulat- 
ing the size of the windows, in which the full time series is 
being decoupled: the transition point is not smeared over 
a wider K t interval, but conversely, it remains strictly po- 
sitioned in K%. Notice however, that the value of K% was 
obtained keeping fixed the geometry of the experiment, 
in particular the radar-target distance: it would be inter- 
esting to vary this distance and monitor the consequent 
changes on the estimated value of K%. An experimental 
study on this specific point is already being planned, and 
results will be reported in a further publication. These 
latter results entail the possibility of establishing a con- 
trol flag in wind intensity, when the device is being op- 
erated in weather conditions which are likely to strongly 
influence the data acquisition process. Indeed, measure- 
ments fluctuating according to skewed distributions could 
in principle induced biased conclusions (e.g. position of 
the mean) without proper information. Moreover, it can 
be speculated that the fitting parameter a returns an 
indirect measure of the correlation degree among air vol- 
ume cells, in which the emitter-detector distance can be 
ideally segmented. When the wind fluctuations are pro- 
nounced, correlations become certainly crucial over the 
finite emitter-detector separation distance and a — > ir/2, 
as seen in 2l|. Qualitative agreement is also found in 



the opposite limit, when the wind turns weak. 

Summing up, in this study it is demonstrated that the 
fluctuations in time of a radar signal propagating through 
the atmosphere, even over short ranges (less than 100 m) 
and moderate wind conditions, are generally distributed 
according to non Gaussian PDFs. Moreover, they gen- 
erally agree with the GG distributions described in ([3]). 
These distributions, have been recently shown to describe 
very well the fluctuations of global quantities in corre- 
lated systems, turbulent systems and many other natural 
phenomena. The conjecture of wind driven fluctuations 
in the radar measurements is further strengthen by the 
observation that the shape parameter a of the GG distri- 
bution fitting the radar data PDFs is directly related to 
the TKE of the crossed medium. All these facts clearly 
support the idea that the non-Gaussian fluctuations of 
radar data are determined by the wind strength. The de- 
veloping of a consistent model for radar signal transmis- 
sion through a windy atmosphere, according to which the 
GG distribution emerges from the physical mechanisms 
of scattering, will be the focus of forthcoming works. 
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FIG. 1: The radar sensor at the outdoor experimental facility. 
The site was equipped with a trihedral corner reflector (high- 
lighted on the right side of the picture) , posed at about 60 m 
from the radar, and an ultrasonic anemometer (highlighted 
on the left). 
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FIG. 2: The PDF of the series A dct (t) (dot line) and the 
best-fit GG distribution (solid line) shown in the upper pan- 
els refer to three different analysis windows. The value of 
the shape parameter a as resulted from the fitting procedure 
is also reported on each panel (notice that G a depends on 
the gamma function, which limits the maximum accessible 
value for the shape parameter to a = 143). The shape of the 
distribution varies from Gaussian, for low TKE values (up- 
per panel, left), to Bramwell-Holdsworth-Pinton distribution 
[111 ], for high TKE values (upper panel, right). The general 
agreement confirms the adequacy of the GG distribution for 
describing the radar data fluctuations. 

Distributions with shape similar to the ones reported in the 
panels are found for TKE values within the intervals high- 
lighted below the panels. This suggests a relationship be- 
tween the shape, and hence the a parameter values, and the 
wind fluctuations intensity. The scatter plot in the main 
panel of the figure reports the inverse of the shape param- 
eter, a -1 , vs. TKE. The direct driving of the turbulent wind 
in the medium toward asymmetric statistics, is stressed by 
the solid line representing a fit based on the sigmoid profile: 
Ai + A2 tanh [A3 (a; — A4)] with A; free parameters. 



